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^ ■ Abstract 

Ch \ Similarities between fragile glasses and spin glasses (SG) suggest the study of 

I ' frustrated spin model to understand the complex dynamics of glasses above 

the glass transition. We consider a frustrated spin model with Ising spins 
and s-state Potts spins both with and without disorder. We study the two 
C^ , models by Monte Carlo simulations in two dimensions. The Potts spins mimic 

orientational degrees of freedom and the coupled frustrated Ising spins take 
F^ I into account for frustrating effects like the geometrical hindrance. We show 

^1 ' that in this model dynamical transitions and crossovers are related to static 

transitions. In particular, when disorder is present, as predicted and veri- 
fied in SG, a dynamical transition between high-temperatures exponential to 
low-temperatures non-exponential correlation functions numerically coincides 
"^ \ with the ordering temperature of the ferromagnetic regions, i.e. the Griffiths 

Q>^ ' temperature Tc{s), while a crossover between power law growth of correlation 

'nI" ■ times to Arrhenius law occurs near a Potts transition at Tp{s) < Tc{s). In the 

f— ^ . model without disorder, where Tc{s) is not defined, both dynamical transition 

and crossover occur at Tp{s). Furthermore the static susceptibility and the 
autocorrelation times of quantities depending on Potts spins diverge at Tp{s). 
C^ ' This is reminescent of recent experimental results on glass- forming liquids. 
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I. INTRODUCTION 



Glassy systems, like supercooled liquids, polymers, granular material, vortex glasses, 

ionic conductors, colloids, plastic crystal and spin glasses, can be defined as materials char- 

^ ■ acterized by a slowing down of one or more degrees of freedom which prevents the system to 

reach the equilibrium as the temperature decreases (Angell 1984, Gotze 1989, Angell 1995, 

Ediger et al. 1996). 

In particular, spin glasses (SG) are dilute alloys of magnetic atoms in low concentration 
in a non-magnetic matrix where the relevant degrees of freedom are the impurities spins 
interacting via random magnetic couplings. The randomness gives rise to frustration, i.e. 
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a competition that prevents the systems to minimize the global free energy (Mezard et al. 
1987). 

One of the characteristic phenomenon of glassy systems is their complex dynamics well 
above the calorimetric glass transition at temperature Tg. For example there is a dynam- 
ical transition at T* > Tg from high-temperatures exponential to low-temperatures non- 
exponential correlation functions. Another feature is the crossover at Ti with Tg < Ti < T* 
between a high-temperature power law growth of correlation times to an Arrhenius be- 
haviour r ~ explAE/ksT), where AE is the activation energy. This phenomenology is well 
seen in experiments (Angell 1984, Angell 1995, Ediger et al. 1996 Andreozzi et al. 1998). 

Here we show the numerical results on frustrated models with and without disorder, 
the Potts SG (PSG) and the Potts fully frustrated (PFF) model respectively, where the 
dynamical anomalies correspond to static transitions, giving some insights on the relation 
between dynamics and statics in glassy systems and clarifying the role played by frustration 
and by disorder. The model could be also useful to understand some recent experimental 
results on finite-temperature diverging static susceptibility in glassy systems (see for example 
Dixon et al. 1990, Menon and Nagel 1995, Leheny and Nagel 1997). 

II. THE PSG AND PFF MODELS 

Let us consider the Hamiltonian 

H = -sjJ2[6^^^^ie,,S,S, + l)-2] (1) 

where to each lattice site is associated an Ising spin Si = ±1 and a s-state Potts spin 
(Tj = 1, . . . , s. The sum is extended over all nearest neighbor (NN) sites, tij = ±1 is 
the sign of interaction and J is the strength of interaction. When the e^j are quenched 
random variables we will refer to the model as the PSG, while when they are distributed in 
a deterministic way to frustrate each lattice cell we will refer to it as the PFF. For s = 1 
the PSG (PFF) model recovers the Ising SG (FF) model. 

The model can be interpreted as a system with s-states orientational degrees of freedom 
frustrated by means of coupled Ising spins which take into account for effects due to geo- 
metrical hindrance (e.g. like in plastic crystals, where the centers of mass of the molecules 
form a regular crystal but the molecules are frustrated respect to the orientational degrees 
of freedom, or in o-terphenyl or in glycerol). The way in which the frustration is introduced 
distinguishes this model by other models, e.g. the Potts glass (Kirkpatrick and Wolynes 
1987, Thirumalai and Kirkpatrick 1998), used to study structural glasses. 

Both PSG and PFF models can be mapped on a Fortuin-Kasteleyn (FK) percolation 
model (Fortuin and Kasteleyn 1972, Coniglio and Klein 1980), defining the FK clusters 
as maximal sets of spins connected by bonds activated between NN spins with probability 
p = 1 — q-'^^J/'^bT ^^^ jg ^YiQ Boltzmann constant) when both the NN Potts and Ising spins 
minimize the coupling energy and using the relation 

Z{e.A= E e-"/'-^ = j:W.,{C) (2) 

{Si,iJi} c 



where Ws{C) = if the cluster configuration C contains any frustrated loop (defined below), 
otherwise Ws{C) = pl'-^l(l — p)l"^l(2s)^'''"^ where N(C) is the number of clusters in the 
configuration C, \C\ is the number of activated bonds and \C\ + \A\ is total number of 
interactions. A frustrated loop is a closed path that has an odd number of antiferromagnetic 
interactions (Coniglio et al. 1991). While the Hamiltonian (0) is defined only for integer 
s > 1, the partition function (§) is meaningful for any values of s (Fierro et al. 1999). 

III. STATICS AND DYNAMICS OF THE MODELS 

The numerical phase diagrams of PSG and PFF model are shown in (Franzese and 
Coniglio 1998, Franzese 1999) and are qualitatively reproduced in figure |I| for the 3D case. 
Two transition temperatures are seen. The lower temperature Tsg{s) (or Tff{s)) corre- 
sponds to a SG (or FF respectively) transition in the universality class of ± J Ising SG (or 
Ising FF) model. The higher Tp[s) corresponds to the percolation transition of FK clusters 
and for any integer s > 1 marks a real thermodynamic transition in the universality class of 
a ferromagnetic s-state Potts model (Wu 1982). At Tp{s) the susceptibility of Potts spins di- 
verges (Franzese and Coniglio 1998, Franzese 1999). In the PSG model, above Tp{s), there is 
the Griffiths temperature Tc{s) defined as the critical temperature of a ferromagnetic model 
with the same number of states of the disordered model (2s for the s-state PSG model). 
This transition is consequence of the presence of ferromagnetic regions, due to disorder, and 
vanishes for vanishing external field (Griffiths 1969). 

From the dynamical point of view we have studied (Franzese and Coniglio 1999, Franzese 
1999) for s = 2 both for PSG and PFF the time dependent nonlinear susceptibility 
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(where A^ is the total number of spins and to is the equilibrium time) that converges asymp- 
totically for t -^ oo to the usual static nonlinear susceptibility. The normalized correlation 
function is 

r- ^ XsGJt) - XsGJt = oo) 

"" XsG{0)-XsG{t = oo) ^ ' 

with Xs'G'(O) = A^. Following (Campbell and Bernardi 1994) the infinite size behavior of /(t) 
has been extrapolated at every t plotting the data for finite linear sizes L = 20, 25, 30, 40 
vs 1/L. 

To test the form of f{t) we fitted the data i) with a simple exponential, finding good fits 
only asymptotically for long time and for high temperatures, ii) with a stretched exponential 
function /o ■ exp[(t/r)^], finding that it fails to fit the data only for short times, ii) with 
the form /o ■ t~^'exp[(t/r)^] suggested by Ogielski (Ogielski 1985), finding that it fits very 
well the data over all the time's range and the temperature's range (Franzese and Coniglio 
1999). 

The parameters /3, r used in the fit ii) and Hi) for xsg are plotted in figures 0, |] for 
PSG and PFF respectively. In the same figures is shown also the integral correlation time 
defined as 



1 t-max 

where / is the generic correlation function. 

In figure ^ for PSG it is possible to distinguish three dynamical regions. A high temper- 
atures region for T > T* c::^ Tc{s) where the correlations functions are exponential {(3 = 1); 
an intermediate region for T* > T > Ti ~ Tp{s) where f3 < 1 and the correlation time's 
growth is a power law forecasted by the Mode Coupling Theory (MCT) (Gotze 1989) or 
mean field (MF) theory (see for example Bouchaud et al. 1998 and references therein); a 
lower temperatures region for T < Ti where activated processes dominate the dynamics (i.e. 
the correlation times follow an Arrhenius law) as in the experimental glasses (Angell 1984, 
Angell 1995, Ediger et al. 1996, Andreozzi et al. 1998). By the way, these results show that 
the prediction T* = Tc given for the SO (Randeria et al. 1986, Cesi et al. 1997) is valid also 
in this model. 

In figure |^ for PFF, where Tc{s) is not defined for the absence of disorder, the intermediate 
region disappears and both the dynamical transition at T* and the crossover at Ti occur at 
a temperature numerically consistent with Tp{s). 

For the PFF model we have calculated also the correlation functions for the Potts order 
parameter M = [s maxj(Mj) — l]/(s — 1) (where Mj is the density of Potts spins in the i-th 
state) defined as 



/M(t) 



{6{t + to)6{to)) 



(6) 



(where 6{t) = M{t) — (M) and to is the equilibration time). At the Potts transition tempera- 
ture Tp{s) the static susceptibility related to M diverges (Franzese 1999) and the parameter 
(3 and r for the fit forms ii) and iii) of Jm becomes less than one and diverging respectively, 
as shown in figure ^. The divergence of r comes out from the observation that t{T, L) grows 
for increasing L at any T and there is a cusp at a temperature consistent to Tp{s) within 
the numerical error. 



IV. CONCLUSIONS 

PSG and PFF model reproduce the experimental phenomenology of fragile glasses above 
the glass transition, recovering the dynamical transition between exponential and non- 
exponential correlation functions and the crossover between the power law growth, pre- 
dicted by MF theories and MCT, and the activated Arrhenius regime for the correlation 
times. In particular, the disordered PSG model shows three different dynamical regimes 
experimentally seen for spin 1/2 probes in glass-forming liquids like o-terphenyl (Andreozzi 
et al. 1998). 

The models show that these dynamical transitions and crossovers can be related to 
static thermodynamic transitions: in the PSG to the Griffiths transition and to the Potts 
transition, while in the PFF model only to the Potts transition. The mapping of these 
models on FK percolation model, show that the Potts transition coincides with a percolation 
transition where frustrated loop are excluded, giving a geometrical interpretation of the 
effects of frustration. 



Furthermore at the Potts transition the susceptibihty associated to the Potts variables 
diverges. This result is reminescent of recent experimental findings (see for example Dixon 
et al. 1990, Menon and Nagel 1995, Leheny and Nagel 1997). In the models presented 
here this divergence is explained by the ordering of the orientational degrees of freedom. If 
such a divergence is present also in structural glasses like plastic crystals is matter of debate 
(Schneider et al. 1999). 

Moreover the PFF model emphasizes that dynamical anomalies are not related to the 
presence of frustration and disorder, but also to the sole frustration and in this case Tp, that 
is the ordering temperature of the orientational degrees of freedom, marks the onset of these 
anomalies. This could be relevant in the study of dynamical behavior of experimental FF 
systems like Josephson junction arrays. 
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FIGURES 

FIG. 1. Qualitatively phase diagram for the PSG model in 3D as function of s. Solid lines mark 
real thermodynamic phase transitions, while dotted line marks the vanishing Griffiths transition. 
Cross marks Potts vanishing transition in the s = 1 (±J Ising SG) case. An analogous phase 
diagram holds for the PFF model with a FF Ising phase at the place of the SG phase and without 
the G phase. 

FIG. 2. PSG for s = 2: (3 and r parameters used to fit the nonlinear susceptibility correlation 
function for L ^ oo (from Ref. [breve] ) . Circles are the parameters for the form in), triangles for 
the form ii); squares are the integral correlation times. Where not shown the errors are smaller 
then symbols size. Arrows show ksTc/J = 3.641 and ksTp/J = 2.925 ± 0.075. 

FIG. 3. PFF with s = 2: As in the previous figure. Arrow shows ksTp/J = 2.76 ± 0.05. 

FIG. 4. PFF with s = 2: as in the previous figure but for the Potts order parameter's 
correlation function Jm- 
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